Magnetic anisotropy switching in (Ga,Mn)As with increasing hole concentration 
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We study a possible mechanism of the switching of the magnetic easy axis as a function of hole 
concentration in (Ga,Mn)As epilayers. In-plane uniaxial magnetic anisotropy along [110] is found 
to exceed intrinsic cubic magnetocrystalline anisotropy above a hole concentration of p — 1.5 x 
10 21 cm -3 at 4 K. This anisotropy switching can be realized by post-growth annealing, and the 
temperature-dependent ac susceptibility is significantly changed with increasing annealing time. On 
the basis of our recent scenario [Phys. Rev. Lett. 94, 147203 (2005); Phys. Rev. B 73, 155204 
(2006)], we deduce that the growth of highly hole-concentrated cluster regions with [110] uniaxial 
anisotropy is likely the predominant cause of the enhancement in [110] uniaxial anisotropy at the 
high hole concentration regime. We can clearly rule out anisotropic lattice strain as a possible origin 
of the switching of the magnetic anisotropy. 
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I. INTRODUCTION 

Ferromagnetism in III-V magnetic alloy semiconduc- 
tors has been studied in many theoretical and experimen- 
tal approaches, Q, 0, 0, |H @| invoked to p — d exchange 
interaction between hole carriers and magnetic impuri- 
ties. One of the most fascinating magnetic features in 
(Ga,Mn)As epilayers is the temperature variation of the 
magnetic anisotropy from (100) to [110] with increasing 
temperature, ja, S 13 where the origin of [110] uniax- 
ial anisotropy remains still unclear. Also, the mechanism 
of the anisotropy switching from (100) to [110] with in- 
creasing temperature had not been elucidated, although 
theoretical works^, have made an attempt to corre- 
late the magnetic anisotropy with the shape of the Fermi 
surface of valence hole subbands. 

In regard to this point, we have recently proposed an- 
other two-phase scenario to explain the magnetic prop- 
erties, where two ferromagnetic phases of [110] uniax- 
ial anisotropic clusters and (100) matrix coexist in a 
(Ga,Mn)As epilaver:[Tol Hll] [110] uniaxial anisotropic 
clusters have relatively high hole concentration in the 
epilayer. Using this model, we can clearly explain the ori- 
gin of the temperature-dependent anisotropy switching: 
with increasing temperature, the ferromagnetic (100) ma- 
trix phase becomes paramagnetic, and, the [110] uni- 
axial anisotropic clusters dominate the magnetism of 
(Ga,Mn)As epilayers near the Curie temperature. [Tol ITl| 
This cluster/matrix scenario is based on the fact that the 
low-temperature grown (Ga,Mn)As epilayers have point 
defects such as Mn interstitials (Mnr), 0, ^| and the 
Mni donors are diffused by post-growth low-temperature 



annealing. [lH IT5I ] Also, the presence of energetically fa- 
vorable Mn clusters, e.g., MnQa-MnQa and MnGa-Mm- 
MriQa, is proposed by theoretical studies [T^.[l7| and di- 
rect observation using scanning tunneling microscopy. 
Recently, Sato et al. theoretically suggested the presence 
of spinodal decomposition phases in magnetic semicon- 
ductors such as (Ga,Mn)As and (Ga,Mn)N, where the 
phases bring about a high Curie temperature. [Tgj 

On the other hand, the effect of hole concentration 
(p) on the magnetic anisotropy has also been exam- 
ined experimentally. For example, it has been reported 
that (100) magnetocrystalline anisotropy can be con- 
trolled by changing p in micro-structured (Ga,Mn)As 
wires^ and modulation-doped (Ga,Mn)As/(Ga,Al)As 
heterostructures . |2 1) indicating that p is a predominant 
parameter in determining the cubic magnetocrystalline 
anisotropy. In terms of in-plane uniaxial anisotropy 
along [110], we have demonstrated that the [110] uni- 
axial anisotropy is enhanced in an annealed (Ga,Mn)As 
epilayer compared with the corresponding as-grown 
epilayer. [2^ Also, Sawicki et al. 23] suggested that the 
in-plane uniaxial anisotropy depends on p using the p — d 
Zener model on the assumption of a small trigonal dis- 
tortion. Since such trigonal distortion is likely to be at- 
tributed to the distribution of Mn ions (Mn clusters), 
this assumption seems to be compatible with our sce- 
nario: the Mn clusters cause inhomogeneity of p in an epi- 
layer and the magnetic anisotropy. To obtain a common 
view about the p dependence of the magnetic anisotropy, 
we hence need to reconsider the magnetic anisotropy 
problem as a function of p on the basis of the cluster 
scenario. [1^ lll| 
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FIG. 1: M ~ H curves of (Ga,Mn)As epilayers at 4 K. Mag- 
netic field is applied parallel to either (100), [110], or [1101. 
Hole concentrations (p) of three samples are (a) 3xl0 20 cm -3 , 
(b) 7xl0 20 cm" 3 , and (c) 1.5xl0 21 cm" 3 . 



In this work, we show a clear manifestation of holc- 
concentration-dependcnt cubic and uniaxial anisotropy 
in (Ga,Mn)As epilayers, and find switching of the mag- 
netic easy axis from (100) to [110] in the high hole con- 
centration regime (p >10 21 cm -3 ) at 4 K. Comparison 
between ac susceptibility data of two separate epilayers 
with p ~10 20 cm -3 and p ~10 21 cm -3 reveals that [110] 
magnetic anisotropy arises from highly hole-concentrated 
(Ga,Mn)As clusters in the (Ga,Mn)As epilayer and such 
clusters are becoming predominant with increasing hole 
concentration. 



II. SAMPLES AND EXPERIMENT 

(Ga,Mn)As epilayers with a thickness of 100 or 200 
nm were grown on a semi-insulating GaAs (001) sub- 
strate using low-temperature molecular beam epitaxy 
at 190 - 240°C. Prior to growing the (Ga,Mn)As, a 
400-nm-thick GaAs buffer layer was grown in an As- 



stabilized condition at 590° C. The Mn contents x in 
the Gai-zMn^As epilayers were estimated to range from 
0.024 to 0.088. After growth, some of the epilayers were 
annealed at 25O°C0, 0, 01 Structural characteriza- 
tions were performed by high resolution x-ray diffrac- 
tion (HRXRD) using a Philips MRD. Magnetization was 
measured with a superconducting quantum interference 
device (SQUID) magnetometer. Ac susceptibility mea- 
surements were carried out using a physical property 
measurement system (PPMS). The hole carrier concen- 
trations (p) which correspond to the ionized Mn accep- 
tor concentration were measured with an electrochemi- 
cal capacitance- voltage method at room temperature. 
The samples we use show typical hole-mediated ferro- 
magnetic properties; Curie temperature (T c ) and hole 
concentration (p) range from 55 K to 135 K and from 
3xl0 20 cm" 3 to 1.5xl0 21 cm" 3 , respectively, being sim- 
ilar to those used in previous work reported by Ku et 
al. 14] 



III. RESULTS AND DISCUSSION 

Figure 1 shows the field-dependent magnetization 
(AI — H curve) of (Ga,Mn)As epilayers with different hole 
concentrations of (a) p — 3xl0 20 cm" 3 , (b) p = 7xl0 20 
cm" 3 , and (c) p = 1.5 x 10 21 cm" 3 at various field orienta- 
tions at 4 K. For sample (a), no differences in the M — H 
curves are seen for field directions along [110] and [110] 
(Fig. 1(a)). In contrast, for samples (b) and (c), the 
differences between [110] and [110] become clearer with 
increasing hole concentration (Figs. 1(b) and (c)), which 
is evidence of [110] uniaxial anisotropy. Also note that 
the remanent magnetization for [110] is larger than those 
for (100) in Fig. 1(c), indicating that the magnetic easy 
axis alters from [100] to a direction near [110] in the high 
hole concentration regime p >10 21 cm" 3 . 

To deduce the magnetic anisotropy constants, we as- 
sume that the M — H curves along the hard axis can be 
modeled on coherent rotation. The magnetostatic en- 
ergy of in-plane magnetized (Ga,Mn) As is generally given 
by E = if u sin 2 (^-45°) + (Jf i/4)sin 2 2^ - MHcos(tp - 
6?), where K u and K\ are the in-plane uniaxial and cubic 
anisotropy constants, respectively, M is the magnetiza- 
tion, H is the applied field strength, ip is the magnetiza- 
tion direction, and 9 is the applied field direction with 
respect to GaAs[100].|l|l|2| We fit a calculated M-H 
curve in the field region from Oe to 2 kOe using the con- 
ditions dE/dtp = 0, d 2 E/dtp 2 > 0, and 9 = 135° ([HO]) 
in this formula. [23 Figure 2 (a) plots K u and K 1 at 4 K 
as a function of p. K\ decreases from 3xl0 4 to 6xl0 3 
(erg/cm 3 ) with increasing p. In contrast, K u gradually 
increases from lxlO 3 to 8xl0 3 (erg/cm 3 ) and becomes 
larger than K\ for a sample with p = 1.5xl0 21 cm" 3 , 
showing significant [110] uniaxial anisotropy. Our exper- 
imental data in Fig. 2 (a) clearly show a crossover of 
the values of K u and K\ in the high hole concentration 
regime p >10 21 cm" 3 . 
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FIG. 2: (a) K u and K\ estimated from M — H curves at 4 K 
as a function of room-temperature hole concentration, (b) K\ 
vs Mn content at 4 K. Different symbols are used for as-grown 
and annealed samples. 
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Although p seems to be the most appropriate pa- 
rameter determining the magnetic anisotropy, another 
cause such as Mn content and lattice strain should be 
checked whether they could represent the experimental 
behavior. Here, Mn content is estimated from HRXRD 
measurements. [28| As shown in Fig. 2 (b), it is clearly 
seen that the plots of K\ vs Mn content are rather scat- 
tered and show no systematic changes as a function of Mn 
content: even samples with the same Mn content show 
different K\ values depending on whether the sample is 
as-grown or annealed. Also, the magnetic anisotropy of 
(Ga,Mn)As was reported to be affected by the shape of 
valence hole subbands and lattice strain could be a pos- 
sible cause of the change in the shape of the valence hole 
subbands. 0, 0| We have confirmed that all the epilayers 
used were grown coherently with tetragonal compressive 
strain and no misfit dislocation, and that no in-plane 
crystal asymmetry between [110] and [110] is induced 
by increasing Mn contents or low temperature anneal- 
ing, ensuring that the influence of the asymmetric lattice 
strain on the change in the values of K u can be ignored. 
Welp et al. also showed that [110] uniaxial anisotropy 
field did not depend on epilayer thickness and the effects 
of (Ga,Mn)As surface and (Ga,Mn)As/GaAs(001) inter- 
face on [110] uniaxial anisotropy are not significant. |29j 
These experimental facts support our description that p 
is most likely the predominant parameter in determining 
the magnetic anisotropy of (Ga,Mn)As rather than other 
parameters. 



FIG. 3: Zero-field cooled and field cooled magnetizations for 
samples with various annealing time ((a) as-grown, (b) 1-hour 
annealed, (c) 16- hour annealed). 

To gain additional insight into the origin of the change 
in the magnetic anisotropy with p, we examine the mag- 
netic features of as-grown and low-temperature annealed 
samples (1-hour annealed and 16-hour annealed samples) 
with a constant Mn content. For the sample with a 
Mn content of 4.4%, the [110] uniaxial anisotropy at 
4 K and the p were enhanced with increasing anneal- 
ing time, being consistent with previous studies. p2l |23| 
Also, no anisotropy switching with increasing tempera- 
ture was observed for 16-hour annealed sample because of 
the marked [110] uniaxial anisotropy. ^l| For these sam- 
ples, we measure the zero-field cooled (ZFC) and field 
cooled (FC) magnetizations. Figure 3 presents the ZFC 
and FC magnetizations of (a) as-grown sample, (b) 1- 
hour annealed sample, and (c) 16-hour annealed sample 
at hdc = 20 Oe along [100]. Clear bifurcations of the 
ZFC and FC magnetizations are seen for all the sam- 
ples. If there are no inhomogeneous magnetic phases 
and the magnetic structure consists of a single-domain 
as pointed out by Wang et at, 30] we cannot expect 
the bifurcation of the ZFC and FC magnetizations for 
the samples, ^l) Thus, these bifurcations strongly sup- 
port the presence of inhomo gen eous magnetic phases in 
(Ga,Mn)As epilayers. [13 [Hill [H EE! 

For our samples, since the single-domain model pre- 
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sented by Wang et aL[30| is not realistic, we hereafter dis- 
cuss the cluster /matrix model in our recent reports: [Tol 
ITlj (Ga,Mn)As epilayers that show two different (100) 
cubic anisotropy and [110] uniaxial anisotropy are com- 
posed of a matrix phase with relatively low hole concen- 
tration and highly-hole concentrated (Ga,Mn)As clusters 
embedded in the matrix. In this model, the distribution 
in p is likely due to the inhomogeneous distribution of Mn 
ions. HH HI E 0E E3 In Fig. 4, we also show the 
temperature-dependent ac susceptibility (X2" —T) curve 
for the same samples as in Fig. 3. An ac magnetic field of 
6 Oe is small enough because features seen in the higher 
temperature regime are compatible with the magnetiza- 
tion data measured in Fig. 3: the high-temperature fea- 
tures in Fig. 4 are seen at around the temperature where 
the magnetization disappears in Fig. 3. As shown in 
our previous study, [10( two peak-to-dip features at ^38 
K and ~57 K in the x^—T curve are seen for the as- 
grown sample, arising from the matrix and the clusters 
with different Curie temperatures . [ToL ITlj Hence, the fea- 
ture seen in the lower temperature regime (^38 K) is as- 
sociated with the Curie temperature of the matrix, (lflj 
For annealed samples, significant changes in the X2" are 
found: the sample subjected to 1-hour annealing shows a 
shift of the low-temperature feature at 38 K towards fur- 
ther lower temperatures, while the high-temperature fea- 
ture at 57 K shifts towards higher temperatures, causing 
multiple-dip feature in the X2 1 ~T curve. After 16-hours 
annealing, the high-temperature feature in the X2" —T 
further changes and shifts towards higher temperature, 
and the low-temperature feature almost disappears. 

We discuss the above systematic changes in the X2" 
— T curves, together with recent studies based on the 
presence of Mn clusters. Mahadevan et al. 16] theo- 
retically suggested that Mn clusters, i.e., MnQa-Mni- 
MnQa complex, are formed favorably under As-rich low- 
temperature growth conditions. Edmonds et oi.[l5| also 
proposed that the dissociation process of Mni-MnQa pairs 
is likely promoted in the region where MnQ a concentra- 
tion is rather high and further low-temperature anneal- 
ing diffuses the dissociated Mni ions. Since the number 
of Mni ions compensating holes decreases with anneal- 
ing in such a region (highly hole-concentrated cluster 
region), the p should increase with annealing, enhanc- 
ing the Curie temperature as seen in the shift of the 
high-temperature feature in the X2" — T curve towards 
higher temperature. It is also likely that the volume of 
the highly hole-concentrated cluster region increases with 
annealing at the same time. On the other hand, some of 
the dissociated Mn ions which diffuse in the annealing 
process may be trapped again in the matrix phase and 
the dissociated Mn ions compensate holes, resulting in 
a decrease in the p in the matrix. Therefore, the low- 
temperature feature in the X2" —T curve is suppressed 
with annealing and simultaneously the contribution of 
the cubic anisotropy which we attribute to the matrix 
phase is also reduced. The other dissociated Mn ions 
may segregate at the surface regionfTEj and have no con- 
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FIG. 4: Temperature dependent nonlinear ac susceptibility 
of the same samples used in dc magnetization measurements 
(Fig. 3) at various ac frequencies. 

tribution to the change in the magnetic properties of the 
matrix phase. As a result, the increase in the p in the 
cluster regions is more significant than the decrease in 
the p in the matrix. These descriptions seem to be con- 
sistent with theoretical and experimental results on the 
Mn distribution in low-tem pera ture grown (Ga,Mn)As 
epilayers. [13, El El El El Hi El El E3 In this con- 
text, the multiple-dip feature seen for 1-hour annealed 
sample is caused by an increase in the p distribution, 
giving rise to the inhomogeneous magnetic properties. 

From these consideration, we conclude that highly 
hole-concentrated (Ga,Mn)As clusters formed during 
the growth are the origin of [110] uniaxial magnetic 
anisotropy on the basis of our recent picture of the clus- 
ter/matrix magnetic structure in (Ga,Mn)As. That is, 
the switching of the magnetic anisotropy as a function 
of p originates from the change in the relative volume of 
highly hole-concentrated (Ga,Mn)As clusters. 

IV. CONCLUSION 

In conclusion, we have observed a systematic change in 
K u and K\ as a function of hole concentration (p) , where 
the contribution of the [110] uniaxial anisotropy becomes 
larger than the cubic magnetocrystalline anisotropy at 
the high hole concentration regime (p ~ 10 21 cm -3 ) at 4 
K. We have revealed that the contribution of [110] uni- 
axial anisotropy is becoming significant with increasing p 
in the highly hole-concentrated (Ga,Mn)As clusters and 
the growth of the clusters. This mechanism is the most 
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probable origin of the switching of magnetic easy axis 
with increasing p. 
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